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Abstract 

The large-scale curvature perturbations induced by spectator anisotropic stresses are 
analyzed across the matter-radiation transition. It is assumed that the anisotropic stress 
is associated with a plasma component whose energy density is subdominant both today 
and prior to photon decoupling. The enforcement of the momentum constraint and the 
interplay with the neutrino anisotropic stress determine the regular initial conditions of the 
Einstein-Boltzmann hierarchy. The Cosmic Microwave Background observables have shapes 
and phases which differ both from the ones of the conventional adiabatic mode as well as 
from their non-adiabatic counterparts. 
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Supplementary anisotropic stresses may arise before matter-radiation equality for diverse 
physical reasons. For instance, on a phenomenological ground, one is often led to consider 
a putative anisotropic stress associated with the dark- matter component [1]. Along a sim- 
ilar line it is also legitimate to consider the possibility that the dark-energy component is 
endowed with an appropriate anisotropic stress [2]. The two possibilities mentioned so far 
are characterized by a common feature: in both cases the anisotropic stress is attributed 
to species which are dominant (or just slightly subdominant) today. This is the case for 
the dark-matter and for the dark-energy components in the ACDM paradigm whose basic 
parameters will be taken to coincide, for the purposes of the present analysis, with the ones 
determined by analyzing the WMAP 5yr data alone 2 [3, 4]. There exist however physical 
situations in which the anisotropic stress does not belong to species which are dominant 
today. Examples in this direction are numerous and include, for instance, spectator scalar 
fields evolving (without being dominant) during the various stages of the life of the Universe 
or even large-scale magnetic fields (see, for instance, [5] and references therein). There is 
an important caveat to the latter statement: unlike large-scale magnetic fields, gravitating 
anisotropic stresses do not couple to the charged species of the plasma. The gravitational and 
electromagnetic effects associated with weakly coupled plasmas can be used to understand 
pre-decoupling magnetism in the presence of the relativistic fluctuations of the geometry 
and this approach leads to the calculation of what are customarily called magnetized Cos- 
mic Microwave Background (CMB) observables 3 . The purpose of the present paper is to 
investigate, in general terms, the interplay between a supplementary anisotropic stress and 
the other sources of anisotropic stress which are present in the system prior to photon de- 
coupling. If the source of large-scale inhomogeneity resides in a spectator anisotropic stress 
it will be interesting to compute the dynamical evolution of the relativistic fluctuations of 
the geometry as well as the temperature and polarization correlations. It will be argued that 
the "stressed" initial conditions for the CMB anisotropies constitute an intermediate case 
between the purely adiabatic and the purely non-adiabatic initial conditions (see, e.g., [9]). 

The time evolution of the system 4 will be parametrized in terms of the scale factor a(r) 
normalized at equality, i.e. a = a/a eq ; within this parametrization the (total) barotropic 

2 For sake of concreteness the values of the cosmological parameters used in the present paper will be 
taken to coincide with the parameters arising from the 5yr best fit to the WMAP data alone, i.e., using stan- 
dard notations, (O b o, ^ c o, ^A, ho, n s , e rc ) = (0.0441, 0.214, 0.742, 0.719, 0.963, 0.087) where e ro denotes the 
optical depth to reionization and n s the spectral index of (adiabatic) curvature perturbations. 

3 The inclusion of an anisotropic stress associated with a large-scale magnetic field is a necessary require- 
ment for the analysis of magnetized CMB anisotropies. It is though not sufficient since the coupling to 
charged particles does lead both to gravitational and electromagnetic effects whose specific account is crucial 
to deduce both the temperature anisotropies and the polarization power spectra. Background magnetic fields 
also affect the photon-electron scattering in a computable way [6]. 

4 Without loss of generality the space-time geometry will be taken to be conformally flat implying that the 
background geometry can be written as = a 2 (r)?7 Ally where r}^ = diag(l, —1, — 1, —1) is the Minkowski 
metric with signature mostly minus. 
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index and the (total) sound speed can be written, respectively, as 
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The exact solution of the Friedmann-Lemaitre equations in a spatially flat background geom- 
etry implies that a(x) = (x 2 + 2x) where x = t/ti and r is the conformal time coordinate. 
The evolution equations for the fluctuations of the geometry and of the sources shall be 
written in Fourier space and in terms of the rescaled wavenumber k = kr\\ if k < 1 the 
given wavelength exceeds the Hubble radius right before the equality time. In the opposite 
case (i.e. k > 1) the given wavelength is smaller than the Hubble radius at the same time 
T\ = (V2 — l)r cq which is roughly half of the equality time when non-relativistic matter and 
radiation contribute equally to the expansion. In the synchronous coordinate system the 
perturbed entry of the geometry reads 5 
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where h(K,a) and a) describe the scalar fluctuations of the metric and enter also the 
contravariant components of the energy-momentum tensor 6 , i.e. 
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The fluctuations of the energy density and of the total anisotropic stress are 7 
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(Pt + Pt)«V t = K 2 [(p u + p u )a u + J^iPa + P*) a * + (Pi + Pl) a x\, 
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(5) 



where J1r(q;) = l/(a+ 1) and J1m(«) = a/(ct+ 1). In Eq. (4) the density contrasts for 
neutrinos, photons, CDM particle and baryons have been introduced (i.e., respectively, 5 V , 
5^, 5 C and <5b); in Eq. (5) <t u (k, a) accounts for the neutrino anisotropic stress while a x (n, a) 
denotes the anisotropic stress of a putative component which is subdominant before equality; 



5 The conventions employed in Refs. [7, 8] differ from the ones employed here (and match with the ones 
of [5] and references therein). The differences result, in a nutshell, from a different signature of the metric, 
from a different naming of the perturbed degrees of freedom and from a slightly different convention in 
the Rayleigh expansion of the brightness perturbations. Bearing in mind Eq. (2) potential ambiguities are 
avoided. 

6 In what follows the arguments of the various functions will be omitted and it will be assume, when not 
otherwise stated, that they do depend both on a and k. 

7 Using standard notations cj c q = /iq^cO, ^bo = h^bo and wmo = /iq^mo- 
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ct x (k, a) will be referred to the photons since photons (together with neutrinos) constitute 
the dominant component of the plasma before equality (i.e. for a < 1). In Eq. (5) a sum 
over other potential components has been added: these components could account for the 
possible stresses arising either in the dark matter or in the dark energy sectors (see, e.g., [1]) 
but will not be explicitly considered hereunder. In the a-parametrization the Hamiltonian 
and the momentum constraints read, respectively, 

dh 
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To write Eq. (7) the following identity has been employed: 
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where 9 t = iK-Vt] moreover, in Eqs. (7)-(8) R u , R 1 = (1—R V ) and Rb(a) denote, respectively, 
the neutrino fraction, the photon fraction and the baryon-to-photon ratio 
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where 3 counts the massless neutrino families, (7/8) stems from the Fermi-Dirac statistics 
and (4/11) 4 / 3 is related to the different kinetic temperature of neutrinos; the numerical 
expression of R\,(a) follows from the choice of parameters listed before in this paper. The 
remaining two equations stemming from the perturbed Einstein equations can be written as: 
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The evolution equations of the neutrinos obey, 

d6 v 26 v . 2dh dB, 
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where Q denotes the combination (/i + 6£) and T v i denotes the £th multipole of the perturbed 
phase-space distribution of the neutrinos. The evolution equations of the dark-matter sector 
obey instead 

db, 6, ldh dd, 6, . 



The system of photons an baryons can be treated, within an excellent approximation, to 
lowest order in the tight-coupling expansion; in the latter approximation the quadrupole of 
the photons vanishes and the resulting equation is given by 

~da a(R h + 1) ~ 8 v / aTT(i2 b + 1) ' 

(i.,. 1 r)h 
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where # 7 b = # 7 — # 7 b represents the common value of the photon and baryon veloci- 
ties. Covariant conservation of the total energy-momentum tensor implies that a x (K,a) = 
3w x Q x (K,a)/4: as well as Q x (k,cx) = [4/ '(3w x )]B(K)(a j 'ai) 1 " 3 ™* . Needless to say that various 
manipulations can be used, in different physical limits, to obtain approximate expressions. 
For instance, by dropping the velocity contribution in Eq. (13) as well as J 7 ^, an approximate 
evolution equation for a v can be derived from Eq. (11), i.e. 

d 2 a u 5a + 4 da v k 2 £ 8 ,„ „ 

ITT + o 7 ^u 7T + T^7 in + . 2 r T7T ^ + ^ - 0. 18 
aa 2 2a(a + l)aa 15(a + 1) 5a 2 (a + 1) 

While Eq. (18) can certainly be used for very small a with fair confidence, it fails as soon as 
either a or k are of order 1. The forthcoming conclusions will not be based on approximate 
equations but on the consistent integration of the whole system in different regimes. 

There are diverse initial conditions which could be imposed to the whole system. They 
are customarily divided into adiabatic and non-adiabatic [9, 10] (see also [7, 8]). In what 
follows we shall be interested in enforcing the adiabaticity condition by demanding that, at 
the onset of the dynamical evolution, the Hamiltonian and the momentum constraints are 
satisfied. The latter considerations imply that, quite generically 8 

3 3 

8 v (k, a ; ) ~ 5 7 (/c, ai) ~ |#b(«, «i) — 4<*c = ~Rj, «i) = 0, Q x (k, a ; ), (19) 

v {k, ai) ~ 7b («, «i) ^ 9 c (k, an) ~ 0, <t x (k, a ; ) = -^^(/c, a ; ), (20) 

where ai denotes the initial integration time when the neutrino anisotropic stress vanishes 
exactly. We are now interested to see what happens to the neutrino anisotropic stress when 
the initial conditions of the system obey Eqs. (19) and (20). The results of the numerical 
integration are illustrated in Fig. 1 for different values of w x and different values of k. The 
initial conditions are set in the limit a a^cc where a^cc = a dec/a e q denotes the value of a 
at photon decoupling (for the best fit to the 5yr WMAP data alone ad ec — 2.92 as indicated 
in Fig. 1 with the vertical line in the plot at the left). If w x > 1/3 the anisotropic stress 
is driven to zero; if k 1, then <j u (k, a) will not oscillate for a > 1 while in the opposite 



8 It should be borne in mind that, in the synchronous gauge, the condition 8 C is enforced not so much 
because of a property of the initial data but rather to fix completely the coordinate system and to avoid the 
occurrence of known spurious (gauge) modes). 
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Figure 1: The relaxation of the neutrino anisotropic stress for different values of w x and k. 

case (i.e. k < 1) the neutrino anisotropic stress will be oscillating for the same range of 
a. If w x = 1/3 the asymptotic value of a u (K,a) will be, approximately, — ^ / R v a x {k, ad ec ) 
for k < 1; in the opposite case (i.e. k > 1) o v will oscillate around the same asymptote 
reached in the n < 1 case. A similar phenomenon has been discussed, in the past, in the 
case of large-scale magnetic fields which affect the anisotropic stress but which also interact 
with the charged particles [11]. Finally notice that, if w x < 1/3, Q x (K,a) grows with a and 
might even get dominant. This situation is more similar to the one treated in [1] and will 
not be specifically addressed here. From the numerical solution of the system in terms of 
h(K,a) and £(k, a) one can also compute other fluctuations with relevant gauge-invariant 
interpretation such as the curvature perturbations on comoving orthogonal hypersurfaces 
(i.e. 7Z), the density contrast on comoving orthogonal hypersurfaces (i.e. e m ) and also ( the 
curvature perturbations on constant density hypersurfaces (see, e.g. [12]) : 

^ . 2a(a + l)d£ K 2 a 2 adQ ^ a + 1 , , 

(3a + 4) dot 6(a + l) 3 da 3a + 4 

Denoting with a\ 1 the value of a at the onset of the numerical integration, £(k, a) oc 
In (a/cti) provided <t v (k, a{) — > since 9 V and 6* 7 b are both proportional to n 2 a. 

It is interesting to assume now that the whole source of large-scale inhomogeneity resides 
in the spectator anisotropic stresses. Could we get reasonable shapes of the CMB observ- 
ables in the absence of the conventional adiabatic mode? The answer is negative and it is 
contained in Fig. 2 where the temperature autocorrelations (i.e. TT power spectra) and 
the temperature-polarization cross-correlations (i.e. TE power spectra) are computed in the 
most favourable case, i.e. w x = 1/3. In both plots of Fig. 2 the full line illustrates the 
WMAP 5yr best fit (in terms of the WMAP data alone and in the case of the conventional 
adiabatic mode) while the dashed and dot-dashed lines refer to the case of the spectator 
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Figure 2: The TT and TE power spectra induced by the spectator anisotropic stress in the 
absence of the conventional adiabatic mode (dashed and dot-dashed lines in both plots). The 
WMAP 5yr best fit is reported, for comparison, with the full lines in both plots. Double 
logarithmic scale is used in the plot at the left while a linear scale is employed in the plot at 
the right. 

stresses with different spectral indices (one of them coinciding with the one of the stan- 
dard adiabatic mode). It has been assumed that the power spectrum of a x is assigned as 
V a (k) = B(k/k p ) rix ~ 1 . The amplitude has been left adjustable 9 and, for the purposes of Fig. 
2, B = 2.41 x 10~ 9 and k p = 0.002 Mpc -1 . The numerical integration has been carried on by 
using, as initial conditions of the Boltzmann solver, the result of the numerical integration 
of the explicit system discussed above. The matching regime between the two regimes coin- 
cides with a ~ 10 -4 when the tight-coupling between baryons and photons is valid. The first 
regime of evolution dictated by the equations derived here avoids a potentially stiff problem 
if a is initially very small. The Boltzmann solvers is a modified version of [5, 13] which, in 
turn, based on [7, 8]. From Fig. 2 it appears that, as expected, the shapes and phases of 
the TT and TE correlations obtained in the case of the adiabatic mode (full lines in both 
plots) differ from the ones induced by the spectator stresses (dashed and dot-dashed lines in 
both plots of Fig. 2). Spectator stresses can be treated and discussed in conjunction with a 
dominant adiabatic mode; absent the adiabatic mode the shapes of the TT correlations have 
intermediate features between the isocurvature humps of the CDM-radiation mode [9]; the 
situation is also different from the case of the magnetized CMB observables when the adia- 
batic mode is absent [5] (see also [13]). There has been recently some confusing statement 

9 By adjustable we mean that it has been taken to match, approximately, with the observed one. However, 
in spite of possible adjustemcnts in the amplitude and well as in the pivot scale, the shapes of the TT and 
TE correlations cannot be made to coincide neither with the ones of the conventional adiabatic mode nor 
with the typical patterns of the isocurvature modes. 
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in the literature. For instance, in [14] (first reference) it is claimed that large-scale magnetic 
fields provide the same shapes for the TT correlations obtainable in the case of an adiabatic 
mode with appropriate amplitude 10 . This statement is in sharp contrast with the present 
investigation. 

In summary, as anisotropic stresses can be attributed to a dominant fraction of the 
present energy (or matter) density [1] , it is also plausible to speculate that there are stresses 
associated with species which are subdominant today. In a rather general framework, by 
enforcing covariant conservation, the initial conditions of the Einstein-Boltzmann hierarchy 
have been numerically investigated using the normalized scale factor as pivot variable at early 
times. The shapes and phases of the temperature an polarization observables are drastically 
different from the ones induced by an adiabatic mode. Spectator anisotropic stresses lead 
then to an interesting kind of initial conditions which can be constrained using the available 
CMB data in analogy with what has been already done in the case of large-scale magnetic 
fields. This analysis is however beyond the scopes of the present paper. 
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